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Abstract 

We propose a method that allows to take into account detector resolution in the partial wave 
analysis event-by-event fit as a special case. Implementation of the method is discussed and 
the applicability of the method is studied for the J/ip —► K* ± K T —> K + K~ir° and J/ij) —> 
A 2 (1430) ± A'T -> K+K-tt° decays. 

1 Introduction 

The most of partial wave analyses (PWA) are performed in the framework of the maximum likeli¬ 
hood method. For a typical set-up the log-likelihood function is s = — y ■ lnP i5 where i runs over 
selected data events and Pi is the probability to observe an event with measured momenta of final 
particles "pi- The latter is given by Pi = €,(7-1 / f (e<j)d<&, where e,; is the selection efficiency for the 
kinematics of event i, <Ji = <r(pi) is the differential cross section (depending on the fitted param¬ 
eters) and the integral gives the overall normalization factor. To take into account the detector 
resolution one has to introduce a convolution with detector response function: ea —> R (g) (ecr). In 
general case such convolution can not be performed within reasonable CPU time. Here we show 
that for a special case discussed below the detector resolution can be taken into account. 

Without loss of generality the method will be described in application to the J/f) —> —> 

K + K~n° decay, where we are interested in measuring the width of K* and assume that J/xfi is 
produced in a e + e~ collider experiment. We also apply the method to J/ip —> AT2(1430) ± A' :f —> 
K + K~n°. 


2 Method 

Typically the cross section is calculated from measured momenta of final state particles. For the 
following we consider the case when these momenta are taken after the kinematic fit, so that the 
total four-momentum is known precisely and is not smeared. 

The method is based on three main assumptions: 

• We assume that the kinematical variables the cross section depends on can be divided into 
two groups such that in the calculation of the convolution of cross section with the detector 
response function, the cross section dependence on the second group of variables can be 
neglected. The qualitative definition will be given later. 

• The mass resolution in the studied kinematic channel is much smaller than the width of the 
studied resonance (in our example <jkk <F k*, where <jk-* stands for the mass resolution in 
the Kir channel). 

• Monte-Carlo simulation of detector resolution and efficiency is consistent with real data pro¬ 
cessing. 
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From general arguments the differential cross section of the J/—>• K + K~n° decay in the center- 
of-mass reference frame depends on four variables {J/tj) in relativistic e + e _ collisions is produced 
with J z = ±1, so we measure 3x3 particle momenta with overall 4-momentum constraint and 
rotation symmetry along the beam axis) denoted by at, k = 1..4. The Breit-Wigner parts of the 
amplitude depend on invariant masses of any two pairs of final particles (the invariant mass of 
the third pair depends on the first two), for computation convenience we took M 2 (K + -k°) and 
M 2 [K~ 7T°). The other two variables are arguments of the angular parts of the amplitude and we 
do not specify them. 

For an event with the measured momenta q , the convolution reads: 

[R <g> ae)(q) = J R{q' - q, q)e(q')cr(a(q , ))dq l « 

~ e(g) J R(q' — q, q)cr(a(q'))dq', 

where R(q' — q , q) is the detector response function (note, due to the kinematic fit not all momenta 
here are independent and that is taken into account in R{q' — q,q)). For the following we will not 
need the explicit form of R(q' — q,q). Then we take integrals and change variables to a: 


R(q' ~ q, q)a(a(q'))dq r = J R(a' — a q ,q)a(a')da' = 

1 8 2 a(a) 


= ( R(Sa , q) (a(a q ) + ^ ^ ^ 6a k + ^ 
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( 1 ) 


Here averages ( 5a k ) and ( 5a k dai) are calculated using R(Sa,q) as weight. The last step is to take 
into account in equation [I] only variables that cross section strongly depends on. Quantitatively we 
assume that 
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where indexes a and /3 run over the first group of variables (M 2 (K + tt°) and M 2 [K 7r 0 )), ( and £ 
run over the second group. Finally one gets the computation formula: 
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Derivatives and averages are calculated numerically, the averages can be taken as constants in the 
fit. Despite of the chosen variables, this expression can be also used for resonances in the K + K~ 
channel. 


3 Notes on implementation 

The decay amplitude in a particular kinematic channel is a product of angular and Breit-Wigner- 
like parts. So, calculating cross section derivatives in the given approximation, one can vary only 
invariant mass squared keeping the angular part constant, which can be easily implemented in the 
existing PWA code. For our tests the first and the second cross section derivatives were calculated 
in the simplest finite-difference scheme (it required calculation of the Breit-Wigner amplitude part 
in six additional points). 

As resolution is much smaller than the width of the studied intermediate resonance the fit can be 
done iteratively: the first fit is performed ignoring detector resolution effects, averages in equation 
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equation [2] for the obtained solution are calculated, than the fit that takes into account the detector 
resolution is performed. For the cases considered below we saw no improvement of the final results 
if an additional iteration was performed. 


4 Method performance: numerical study 

In our numerical study we assume a typical set-up of a e + e _ collider experiment (e.g. see HD- We 
model the detector response by smearing generated particle momenta. To get an estimate of the 
method sensitivity to particular shape of the detector response, we consider two detector response 
models (referred in the following as “model 1” and “model 2”). In the first model we add Gaussian 
fluctuations with zero mean to track helix parameters (re, tan A, <j>o). For photons the energy and 
the direction are smeared ( e ph , 9 ph , tp ph ) in the same way. To study the method applicability for 
different K ± n° invariant mass resolution we use set generated samples with variances proportional 
to 


(J k /k = 5 x 10 3 , 
^tan A = b X 10 , 

(T Vo = 2.5 X 10~ 3 


for kaon tracks and 


a eP h / e ph = 2.5 x 10” 3 , 
(70ph = 1 X 10 , 

(7 ^pph =i lx 10 


for photons. As we mentioned above, the input for PWA are particle four-momenta after the 
kinematic fit (here the fit is additionally constrained for the 7r° mass). 

In the second model we vary Cartesian momentum projection of A' + , K ~, 7r° with standard 
deviations proportional to 


a n K = a„K = a. n K = 3 MeV, 

r X “y y z 

(7 _o = a _o = a n o = 5 MeV. 

Pi Pi Pi 


Reactions J/ij) —> K*(892) ± K :f —» K + K~tt° and J/ij) —> /V 2 (1430) =l= —> K + K~n° are 

modeled by weighting a phase space distributed sample. For each set of smearing parameters 
the fitting procedure is repeated several times with different generated phase space samples of 
approximately 0.9 x 10 6 events. 

The decay is parametrized in the covariant tensor formalism framework [2] , for the Breit-Wigner 
part we used 


1 

M 2 — s - iMTj(s) ’ 

P(s) r 

p(M 2 j) ’ 

2k k 2J 

yfs F 2 (fc 2 , Tj, J) ‘ 

Here s and k are the invariant mass squared and relative momentum of the resonance daughter 
particles, J is the spin of the resonance; M, F and r are mass, width and Blatt-Weisskopf radius 
correspondingly; functions F(k 2 , r 2 , J) are Blatt-Weisskopf form factors, which can be found in j2|. 
This later is the most suitable parametrization for I\*( 892) and can be also applied for 7\2(1430). 
In our study we use PDG averages [3] for the mass and the width of A"* (892) and A' 2 (1430) and 
fix their Blatt-Weisskopf radii to 0.5 fm. Dalitz plots for the “generated data samples” are shown 
in figure [T] 

Averages (5a a ) and ( 5a a 5ap) are determined iteratively as explained above. 


A bw = 

r As) = 

p{s) = 
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The comparison of fit results that ignore and take into account the correction for the detector 
resolution are shown in figure [2] Deviation of mass, width, Blatt-Weisskopf radius and the fraction 
of corrected width bias are given as a function of the £ = ResMunl T ratio, where RcsMktt is the 
variance of the Mk-k- The invariant mass variance is taken in the resonance band region (i.e. for 
K*{ 892) it is calculated in the band 0.8 GeV < M K + 7T o < 1.0 GeV and M(K~n°) > 1.0 GeV). 
The results are provided for two detector response models. 

5 Discussion 

Firstly, we see that the used method allows to compensate approximately 80-90% of the resonance 
width and Blatt-Weisskopf radius biases due to the detector resolution for £ < 0.2. The method 
applicability dramatically decreases at higher £ values. 

Secondly, due to made approximations the applied method introduces bias to the fitted mass 
of the resonance. In the case of K*( 892) it is almost independent of the detector response model 
and equals to 0.2 MeV for £ = 0.1. If we increase the mass of the resonance this bias decreases 
and finally changes sign. In the case of I\ 2(1430) taking into account detector resolution can 
improve or worsen fit results depending on the detector response model. However, sometimes event 
reconstruction itself can bias measured mass of two particles (for example due to non-Gaussian 
detector response to photons). In the proposed method such effects are taken into account in the 
linear term in equation [2] and possibly can improve overall mass measurement. 

Thirdly, we also apply this method to resonances in the K + K~ kinematic channel and find 
results similar to the presented. 

We found the fit results not essentially dependent on used detector response model, but some 
difference tells us that MC study is needed in any practical application of the method. 

The computation time increased approximately proportional to the number of additional points 
used to calculate the cross section derivatives. 
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Figure 1: a)-b) Dalitz plots for generated J/ip 

/^(IdSO)^A'T —> K + K~tt° samples correspondingly. 


K*(892) ± K^ K + K~ir° and J/if> 


6 Conclusion 

For the first time we propose a practically acceptable method, which allows in special cases to 
take into account the detector resolution in the event-by-event partial wave analysis fit. The 
method reduces the computation of convolution of a process cross section and a detector response 
function to calculating cross section derivatives and tabulating variances of essential cross section 
arguments. We demonstrate the method performance and applicability in the set-up of a typical 
e + e“ experiment for two toy detector response models and two reactions: J/ip —> K*(892) ± K^ : —► 
K + K~n° and J/tp ->• JT 2 (1430) ± it' =F -H- K+K~t: 0 . 
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Figure 2: A comparison of fit results when detector resolution is taken into account (solid markers) 
and ignored (open markers) for K*{ 892). Results for the first and the second (see in the text) detector 
response models are shown in red and black correspondingly. Figures a)-c) show fitted width, mass and 
Blatt-Weisskopf radius, figure d) shows the fraction of width bias due to detector response compensated 
by the used method. The shown results for the model 2 are limited by £ < 0.23 as for higher £ values 
using of the quadratic approximation results in non-positive differential cross section for some “data 
points”. 
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Figure 3: A comparison of fit results when detector resolution is taken into account (solid markers) 
and ignored (open markers) for i\2(1430). Results for the first and the second (see in the text) 
detector response models are shown in red and black correspondingly. Figures a)-c) show fitted width, 
mass and Blatt-Weisskopf radius, figure d) shows the fraction of width bias due to detector response 
compensated by the used method. 
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